Metallicity Dependence of the Cepheid Calibration 
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The metallicity dependence of the Cepheid period-luminosity (PL) relation has been a 
controversial issue. The more commonly accepted view is that the PL relation is insensitive 



to metal abundance, while the period-luminosity-colour relation is fairly sensitive, as indi- 
cated by stellar pulsation models 1-4 (see also refs. 5,6). The direct observational evidence 
is scanty, however. Freedman and Madore 7 have studied Cepheids in three regions having 
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different metallicity in M31 and concluded no statistically significant residual differences 
in the distance moduli. On the other hand, Gould 8 has claimed that the data of Freedman 
and Madore are rather indicating a large metallicity dependence A///A[Fe/H] = 0.88±0.16, 
where (x is the distance modulus derived from the Cepheid PL relation. A moderate metal- 
licity dependence A/i/A[Fe/H] = 0.44 was also concluded in a recent report from EROS 



photometry for the Magellanic Clouds 9 . 

Recently an accurate measurement of metallicity has been carried out for 23 Galactic 
Cepheids by Fry and Carney 10 (hereafter FC). This work has given us an opportunity to 
examine the metallicity dependence of the Cepheid PL relation and its calibrations. The 
most commonly used calibration of the Galactic Cepheids is probably that of Feast and 
Walker 11 based on V band photometry and the distance based on the ZAMS fitting to 
the clusters that contain the Cepheids. Laney and Stobie 12 (hereafter LS) have slightly 
updated the Feast- Walker distance and also provided a calibration of the Cepheid PL 
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relation for the near infrared JHK colour bands 13 in addition to that for the V band 
(their result for V band hardly differs from that of Feast and Walker). The work of LS 
may perhaps be taken as the most accurate calibration of the Cepheid PL relation to date, 
hence giving the most accurate distance to the LMC with the Cepheid method. 

The FC sample contains 12 Cepheids in common with the LS sample: EV Set, V340 
Nor, WZ Sgr, V Cen, SW Vel, CV Mon, S Nor, U Sgr, TW Nor, SV Vul, QZ Nor, and 
T Mon in the order of increasing metallicity. We show in Fig. 1 an uncritical display of 
the residual of the fit to the Cepheid PL relation as a function of metallicity [Fe/H]. The 
residual is defined by 

6M x = M x -(A x \ogP + <t> XPL ), X = V,J,H,K (1) 

with the coefficients given in LS: (A, cfy = (-2.874, -1.197) for V, (-3.306, -1.971) for J, 
(-3.421, -2.243) for H and (-3.443, -2.297) for K (we take the first set of their parame- 
ters, but other choices hardly modify the plot of Fig. 1) . Virtually the same plot is obtained 
with the data and the fit given in Feast and Walker for the V band. The error of the metal- 
licity measurement of FC is ±(0.02 - 0.03) dex except for SV Vul (±0.04) and CV Mon 
(±0.06). This figure shows a conspicuous correlation of M x with metallicity. The surprising 
fact is that this metallicity dependence is disturbingly large 5M x /5[Fe/H] = —(1.3 — 1.4), 
almost independent of the colour band from V to K. We note that one highly deviated 
point (TW Nor) suffers from a very large reddening correction Eb-v (OB) = 1.34 and My 
may not be accurately determined. The leftmost point EV Set and QZ Nor (the second 
point from the right) are suspected to be overtone pulsation (FC). If we remove these three 
Cepheids (this means we remove all highly scattered points in Fig. 1), the metallicity gra- 
dient is 5M x /5[Fe/ti\ = —(1.7 — 2.1). In order to underscore the colour independent nature 
of the offset we plotted SM X (A = V, J, H) versus 5M K in Fig. 2. 

Now the question arises as to whether this large metallicity dependence represents that 
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of the intrinsic Cepheid PL relation or it arises from some other artefacts of procedures 
used in estimating M\. Let us remember that the absolute magnitude is calculated with 

Vo = fi v -R(OB)E B _ V (OB) (2) 

and 

M x = m x - {A x /A v )R(Cepheid)E (Cepheid) - ^ (3), 

where the V band modulus \iy is obtained by ZAMS fitting assuming the cluster metallicity 
[Fe/H]=0 (see e.g., refs. 11, 14). The ZAMS track with lower metallicity generally yields 
fainter magnitude for given colour (e.g., ref. 15) and hence this is just opposite to the 
trend observed in Fig. 1: a metallicity correction for the ZAMS fitting would steepen the 
slope. It is unlikely that the metallicity dependence arises from the absorption correction 
for Cepheids in eq. (3), because the correlation seen in Fig. 1 is wavelength independent, 
whereas reddening is not. Then we are left with two possibilities: 

(i) the Cepheid PL relation indeed has a large metallicity dependence, or 

(ii) the absorption correction for /j, has a large [Fe/H] dependence. 

Case (i), if it does not looks very likely, does not conflict with any other observations, 
when the uncertainties in separating absorption corrections from the metallicity effect and 
also of the metallicity measurements are considered. The indicated metallicity dependence, 
however, is by more than a factor 10 larger than the prediction of the theoretical models 
of the Cepheid pulsation 1-4 . The colour-band independence of the offset at first glance 
also looks a little unusual from a theoretical view point. This is not impossible, however, 
if the coefficient of the pulsation equation (P as a function of mass and radius) is strongly 
affected by metal abundance, in contrast to what stellar model calculations suggest. 1 ' 3 We 
remark that the sign of the period shift is consistent with the Oosterhoff period shift for 
RR Lyr, 16 ' 17 for which stellar models also predict the metallicity dependence by a factor 
of > 10 smaller than observed. 18 If case (i) is true, an on-going Hubble Space Telescope 
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project for studying the metallicity dependence for Cepheids in M101 9 should detect the 
effect; the raw distance moduli derived from inner and outer regions may differ by 1 mag 
after absorption correction. 

In order to examine the model prediction for the metallicity dependence of the period- 
colour (PC) relation, we show in Fig. 3 the residual of the fit of [B— V)q obtained by LS as a 
function of [Fe/H], and compare it with a prediction of Stothers. 3 Unfortunately, the scatter 
of data is too large that one can hardly conclude whether the residual shows the metallicity 
dependence. The only conclusion one can draw from the figure is that the metallicity effect 
on the PC relation, if any, is not too much larger than the theoretical model predicts. 
However, this does not preclude the possibility that the pulsation equation is strongly 
affected by metal abundance: the contribution of the metal-abundance dependent term to 
the zero point of the PC relation is one order of magnitude smaller than that to the zero 
point of the PL relation (see ref. 3). Namely, this test does not serve for us to choose 
among the two possibilities. 

The colour independent nature may also be taken to be consistent with case (ii). 
Namely, either selective extinction i? B _y(OB), which is usually estimated from cluster 
reddenings or space reddenings of stars close to Cepheid, or R factor, which is set to be 
a constant, receives a large abundance effect, although, this differs from our conventional 
belief that the extinction correction for B stars does not depend so strongly on metallicity, 
for line blanketing effect is smaller for high temperature stars. Also somewhat strange with 
this case is that the extinction correction in (1) does not tend to cancel against that in (2) 
for the V band. We have examined possible correlation between E(OB), i?(Cepheid) or 
t^y — Ho an d [Fe/H], but have not found any apparent correlations between those quantities. 

One may suspect a systematic error in FC in estimating metallicity. One of the 
suspects may be in their estimate of temperature. We have checked that their values of 
[Fe/H] do not show any significant correlation with colour (B — V)o, indicating no obvious 
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errors in the temperature estimate that correlates with metallicity. Indirect, but more 
significant is Fig. 2 which demonstrates a strong correlation among residuals of fitting in 
different colour bands, pointing towards a systematic problem in the Cepheid work itself, 
rather than in the estimate of metallicity of FC. (We remark that similar correlations in 
the residuals are noted in refs. 8 and 9.) Namely, possible systematic error of FC, if any, 
does not account for the entire problem posed here. 

Neither of the two cases discussed above looks very likely, but certainly not impossible. 
If (i) is correct, the period- luminosity relation for the Cepheid is represented better with 
the addition of a metallicity term, for example for the V and K bands, such as 

My = -1.256(±0.046) - 2.874 log P - 2. 15(±0.44) [Fe/H] (a = 0.111) (4a) 

M K = -2.385(±0.039) - 3.443 log P - 1.72(±0.38) [Fe/H] (a = 0.095) (46) 

as obtained by fitting our 9 Cepheids (the slope is fixed to be that from LMC; EV Set, QZ 
Nor and TW Nor are excluded) . These are compared with 



without the metallicity term (these equations agree with what are given in LS within the 
error). The equation with the metallicity term yields the LMC distance 0.4 mag smaller 
than estimated in LS for average metal abundance of the LMC, [Fe/H] =-0.3. If (ii) is 
correct the "true" PL relation after correcting for the metal abundance effect on no is 
given by the equation with the third term absorbed into eq. (4) . A comparison of (4) with 
(5) indicate that the LMC distance may be w+0.1 mag larger. This shows that one cannot 
determine the distance accurately from the Cepheid PL relation unless the metallicity 
effect is well controlled. 



My = -1.144(±0.077) - 2.874 log P (a = 0.232) 



(5a) 



M K = -2.295(±0.063) - 3.443 log P (a = 0.189) 



(56) 
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A positive aspect, on the other hand, is that the intrinsic scatter of the Cepheid PL 
relation should substantially be smaller than is believed now once the metallicity correction 
is made. After the metallicity correction the dispersion of the PL relation becomes 0.08- 
0.11 mag depending on the colour band, which is about a half the scatter obtained by 
LS. 

The present authors were unable to penetrate into details of the numbers obtained by 
LS beyond the point which is explicit in their paper. It is, however, obvious that the best 
present data used for the Cepheid calibration, which is a cornerstone of the extragalactic 
distance scale, contains an unusually large metallicity dependence that has been ignored 
in the literature. This would bring a significant uncertainty into the LMC distance; the 
error of the LMC distance from the Cepheid PL relation remains uncertain as large as -0.4 
to +0.1 mag, till this problem is solved. 
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Figure captions 

Fig. 1 Residual of the fit to the Cepheid PL relation as a function of metal abundance [Fe/H]. 
Here the residual is defined by 5M X = M\- (A\ logP + </>apl) (A = V, J, H, K) with 
the coefficients A and 4>\pl) taken from LS together with the data of M\. 

Fig. 2 Correlation of the residuals for different colour bands. 8M V , SMj and 5M H are plotted 
against 8M K . 

Fig. 3 Residual of the fit to the Cepheid period-colour (B — V) relation as a function of 
metal abundance. The residual is defined by 8{B — V)q = (B — V) — (ClogP + D). 
Both coefficients and data are taken from LS. The curve shows S(B — V) ~ 8ASZ as 
predicted by Stothers. 3 
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